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INTERPLAY OF BOLTZMANN EQUATION AND CONTINUITY 
EQUATION FOR ACCELERATED ELECTRONS IN SOLAR FLARES 

A. CODISPOTit AND N. PINAMONTitt 


Abstract. During solar flares a large amount of electrons are accelerated within the plasma 
present in the solar atmosphere. Accurate measurements of the motion of these electrons start 
becoming available from the analysis of hard X-ray imaging-spectroscopy observations. In this paper, 
we discuss the linearized perturbations of the Boltzmann kinetic equation describing an ensemble of 
electrons accelerated by the energy release occurring during solar flares. Either in the limit of high 
energy or at vanishing background temperature such an equation reduces to a continuity equation 
equipped with an extra force of stochastic nature. This stochastic force is actually described by 
the well known energy loss rate due to Coulomb collision with ambient particles, but, in order to 
match the collision kernel in the linearized Boltzmann equation it needs to be treated in a very 
specific manner. In the second part of the paper the derived continuity equation is solved with some 
hyperbolic techniques, and the obtained solution is written in a form suitable to be compared with 
data gathered by hard X-ray imaging—spectroscopy telescopes. Finally, a first validation of the model 
with NASA Reuven Ramaty High Energy Solar Spectroscopic Imager spectrometer measurements is 
provided. 

Key words. Solar flares, Boltzmann equation, Landau equation, continuity equation, electron 
transport. 
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1. Introduction. As seen from the Earth the Sun appears quiet and stable. 
However, complicated phenomena can be observed, looking closer at its surface. 
Among all perturbations occurring in the solar atmosphere, solar flares are the most 
energetic ones. Actually, during solar flares, huge quantities of energy are released. 
In the case of arc-shaped flares, according to the standard picture [8, 17, 38], a large 
number of electrons are first accelerated above the loop top and then injected into the 
loop, finally they move down towards the chromosphere following the magnetic field 
lines. During this motion, electrons lose energy because of various physical processes. 
The most efficient energy loss mechanisms are Coulomb collisions with ambient elec¬ 
trons present in the solar atmosphere [9, 25, 15]. However, these kind of interactions 
are not perfectly elastic and, in every collision, a small, almost negligible, part of 
energy is radiated away due to Bremsstrahlung, mainly in the Hard X-ray domain 
[30, 23]. Photons emitted by means of this mechanism are visible also far away from 
the solar surface. Thanks to this effect we are able in principle to indirectly detect 
electron densities by analyzing hard X-ray photons emitted during flares. 

In the X-regime, solar flares are observed by satellite imaging-spectrometers such 
as, YOHKOH [1], RHESSI [26], which is still in operation, and STIX in the future 
Solar Orbiter mission [5]. Since these instruments are capable of collecting data both 
spatially and spectrally resolved, precise images of photon fluxes at various energies 
can be constructed. In particular, in [32], Piana et al showed how to construct maps 
of electron fluxes E(x, E) averaged over the line of sight and thus information about 
electrons (electron maps) can be obtained from satellite spectrometer data. Nowadays, 
these kind of measurements are becoming very precise and permit to test refined 
models of electron transport mechanisms. 
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Unfortunately, the information about the direction of electrons velocities is lost 
in the detection procedure, because only their energy can be measured. We shall 
see that this issue slightly modifies the validation process of theoretical models by 
comparison with the measured data (i.e., the images provided by these satellites). 

From the theoretical point of view, it is believed that during solar flares electrons 
are accelerated by magnetic reconnection. However, the actual mechanism governing 
this acceleration is still under investigation [4, 10, 31, 42, 6, 20]. For this reason several 
statistical mechanisms of electron acceleration and various theoretical approaches have 
been developed in order to model the physics of solar flares. 

Usually, since electrons are indistinguishable particles, their ensemble is described 
by a density function /(t,x, v) which expresses the number of particles with velocity 
V at the time t in the position x. The integral of / over all possible velocities gives 
the standard particle density n(<,x). The evolution of the density function is given 
by various equations of plasma physics adapted to the kind of plasma present in the 
solar atmosphere. As an example, we recall the analysis presented in [21, 34, 35], 
where perturbations of charged particles in a plasma are shown to satisfy a Fokker 
Planck equation. Recently several theoretical diffusion models were considered and 
adapted to solar flares e.g. [31, 6]. 

Among all of the various possibilities, the continuity equation adapted to the 
phase space was employed to describe the electron motion along magnetic field lines 
in solar flares. We recall that the standard continuity equation has the following form 

(1-1) ^ = |^+vVx/+--Vv/ = ^ 

where a source term S is considered to model injection of electrons in the flaring region 
due to magnetic reconnection. Furthermore, F is the total external force acting on 
electrons and m is the electron mass. Notice in particular that whenever the force F 
is positional the previous equation is equivalent to 

(1.2) ||- + Vx-(v/) + Vv- =5. 

However, in order to take into account the statistical energy loss due to Coulomb 
collisions, the previous two equations are modified by adding the average value of 
a stochastic force F^ to the total external force F. In particular F^ is chosen to 
match the energy loss rate (VxA) computed as the mean value of the energy loss in 
collisions with ambient particles, as e.g. [15, 36]. With this modification the continuity 
equation (1.1) takes the form 

o^ ^ F VT ^ , C^xA) ^ ^ ^ 

(1.3) -^Tv-Vx/H-• Vv/H -• VvJ — o . 

at mm 

A slightly different equation, valid for accelerated electrons propagating in a solar 
plasma, has been derived by Emslie, Barrett and Brown [16] as a consequence of 
conservation of the electron flux along every energy channel. It has been widely used 
in order to understand some details about solar flares e.g. [8, 9, 37]. Furthermore, 
it has recently been employed for the validation of various statistical energy loss rate 
models [39, 12]. 

However, in general the stochastic term Fg = {Vy^E) is not positional, it explicitly 
depends on the velocities, hence equation 
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(1.3) is not anymore equivalent to 


(1.4) 


^ + Vx • (v/) + Vv 



= S 


contrary to the case of purely positional external forces. In this paper, we show 
that only continuity equation in the form (1.3) matches the asymptotic form of the 
Boltzmann equation which is a non-linear kinetic equation that takes into account in¬ 
teraction/collisions of particles by means of a non linear collision kernel. In particular, 
we shall obtain the stochastic contribution in the form 


Lc(/)(t,x,v) = . Vv/ 

m 


(1.5) 


as an asymptotic approximation of the Boltzmann collision kernel. Furthermore, since 
the physical system we have in mind is open, the equation from which we start is the 
Boltzmann equation equipped with some appropriate given source term S', which is 
related to S in (1.1), 


df F 

7^+v-Vx/+--Vv/ = g(/,/) + 5' 
ot m 


( 1 . 6 ) 


where Q{f, f) is the non linear Boltzmann collision kernel for Coulomb collisions and 
where, as before, F are the external classical forces. In this paper we treat Q{f, /) 
in the Landau approximation [13], namely in the limit of grazing collisions, or in 
the limit of large values for the Coulomb logarithm. Furthermore, we consider linear 
perturbations Sf over solutions / of the Boltzmann equation. In the asymptotic 
regime, i.e. for large velocities jvj or for small temperatures in the case of Maxwell 
background, the perturbation Sf satisfies an equation of the form (1.3) provided the 
stochastic force is taken into account with a term of the form • Vvf, which is 

slightly different from that considered in [39, 12]. 

After having obtained the continuity equation (1.3) equipped with the stochastic 
force used to model Coulomb collisions, we shall discuss how its outcome can be com¬ 
pared with measurements. As a first step, we have to adapt the obtained continuity 
equation for electron densities resolved in position and velocities, to describe electron 
fluxes resolved in position and energy. In order to accomplish this task, an averaging 
procedure over the velocity directions needs to be employed. Operating in this way 
we obtain a new continuity equation and the study of the differences between this new 
model and the original one in [16] is at the core of this paper. We will address such 
study from both a theoretical viewpoint and by means of a computational validation 
based on empirical data measured by RHESSI. The paper is organized as follows. In 
the next section we shall recall some standard facts about the Landau approximation 
of the Boltzmann collision kernel for Coulomb interactions and its linearization. In 
the third section we shall analyze the asymptotic form of the linearized Boltzmann 
equation and we establish the connection to the continuity equation employed in the 
literature. In the fourth section we shall discuss the solution of the continuity equation 
in the case of isotropic perturbations. We shall in particular highlight the hyperbolic 
character of the obtained equation. In the fifth section we shall compare the obtained 
continuity equation for electron fluxes with the one presented in [16] also discussing 
their validation with measured data. Hence, we shall show how the proposed modi¬ 
fied equation is consistent with experimental data provided by imaging spectroscopy. 
Finally, some comments are drawn in the last section. 
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2. Landau approximation of the Boltzamnn collision kernel and its lin¬ 
earization. The collision kernel present in the Boltzmann equation (1.6) takes into 
account elastic collisions among particles distributed according to some density func¬ 
tion /(t, X, v). Furthermore, since collisions are supposed to be elastic, the momentum 
and the kinetic energy are conserved in every collision. It means that, if we indicate 
by V and Vi the interacting particles velocities before the collision, and by v' and v( 
the interacting particles velocities after the collision, the following relations hold 

V -k Vl = v' -f v'l, |vp -k |vip = |v'p -k Iv^p. 


From these equations, indicating by cr = (sin0cos^, sin0sin(^,cos0) the deviation 
angle, we can write explicitly v' and v( as 


V -I- Vi 


V- Vi 


, V-Vi |V-Vi| 

Vi — -—- — -—- G. 


The collision kernel Q{f, /), which accounts the interactions between colliding parti¬ 
cles, takes the form [18] 


(2.1) Q(/,/)(t,x,v) = / [ B(|v-Vi|,cr)(/7{-//i)d^crd^Vi 

where, in order to shorten the notation, it is understood that /' = /(v'), /i = /(vi) 
and /{ = /(v(). Furthermore, in the case of Coulomb collisions, the integral kernel 
B is related to the Rutherford cross section 


( 2 . 2 ) 


R(|v - vi|,ct) = 


m2 |v - vi|3 sin"^ | 


o>eo 


where Oq is a cutoff representing the minimum scattering angle, which is related to 
the Debye screening length, and to the Coulomb logarithm A 


A = — log sin - 



The typical values of the Coulomb logarithm for the plasma in the solar atmosphere 
are sufficiently large to justify the assumption of grazing collisions, namely to consider 
only small scattering angles [2]. As a rather direct consequence, we obtain an approx¬ 
imate expansion Ql of Q through a first order expansion in the scattering angle and 
the evaluation of the integral over all possible scattering directions. In this way the 
Landau approximation of the Boltzmann collision kernel is obtained. 

The Landau equation was derived for the first time by Lev Davidovic Landau [24]. 
The limits in which this equation represents a good approximation of the Boltzmann 
equation and the relations between the solutions of the two equations have been widely 
analyzed [13, 14, 7, 27]. In particular, the solutions of the Boltzmann equation are 
well approximated by solutions of the Landau equation if the density function / is 
spatially homogenous [3, 19, 41]. The analogous result in the non homogeneous case 
has been rather recently obtained in [2]. 

The Landau approximation of the Boltzmann collision kernel, which is also called 
Fokker Planck collision kernel [13], has the form 

(2.3) Q^(/,/)(i,x,v) = ^^V- / A(/iVv/-/VvJi)d3vi 

dR3 
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where A is a linear operator from to whose corresponding square matrix has 
the following entries 


(2.4) 


1 


0 - 2-7 - 


4|v - vi| 


Sij 


(v- Vl),(v- Vl)j 
|v - Vi|2 


Suppose / is a solution of the Landau equation 

and consider a small perturbation f + Sf around that solution. 

Since the perturbation Sf is regarded as small with respect to /, we simply de¬ 
duce the equation for Sf by linearization of the Landau equation for /. In particular, 
we neglect terms of higher order in Sf inside the Landau collision kernel, which cor¬ 
responds to the assumption that Q{Sf, Sf) 0. Allowing for a source term S' in 
order to take into account that the system we are considering is open, we find that 
the perturbation Sf is such that 

(2.5) ^ = LBiSf)(t, X, v) -b SBiSf)(t, X, v) -b S', 

at 

where, for fixed /, Ls(5/)(i, x, v) is an operator which acts locally in v on 5/, i.e. 

(2.6) LB(<5/)(t,x,v) = 5^V- / A(/iVv5/-<5/VvJi)dV, 

™ Jm.3 

while Sb is an operator in integral form on Sf 

5B((5/)(f,x,v) = ^^V- [ A((5/iVv/-/Vv,^/i)d3vi. 

TO Jr3 


3. Reduction to the continuity equation. In this section we shall see how the 
linearized Landau approximation of the Boltzmann equation reduces to the continuity 
equation. This result is obtained for large values of |v|, whenever the background 
solution / G C°“(R^) decays rapidly in the velocity domain, namely whenever / at 
fixed time and position is a Schwartz function /(t,x, •) G §(R^). The same outcome 
holds in the case of thermal backgrounds described by Maxwell distributions at small 
temperatures. 

3.1. Asymptotic expansion of the linearized Landau collision kernel. 

We are interested in studying the asymptotic behavior of Sb and Lb for large |v|, 
whenever the background solution / is of Schwartz type in the velocity domain. This 
is the case, for example, whenever / is a Maxwell distribution or a shifted Maxwell 
distribution. 

Let us start discussing Sb in (2.5). Since Sb is a local operator in /, the asymp¬ 
totic behavior of Sb for large values of |v| is governed by /, whenever <5/ is sufficiently 
regular as is assumed through this paper. Hence, since / G S(H), Sb decays faster 
than any inverse power of |v| for large |v|. 

Let’s focus now on the local term LB(t,x,v) in (2.5), where we expect to find 
the dominant contribution for large values of |v|. Integration over all possible target 
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velocities vi, recalling the form of the operator A given in (2.4), provides 


LB(5/)(t,x,v) = - 

m ^ 


- 27r6ff - 


ASf 


h 


/r3 |v- vil 




(3.1) 


9ijSf f (v - vi)*(v - vi)-^ 


/R3 


V- Vi 


fid^vi 


In order to discuss the asymptotic form of all these contributions, we need the following 
theorem 

Theorem 3.1. Let f be a function in S(]R^) and n := fid^vi. The following 
relations hold 


1 . 


2 . 


3. 


f 7 -^—= A+^i(v), 

/k3|v-Vi| |v| 


r (v - vi)^(v - vi)-^ 

/r3 |v-Vi|3 

ft/i ^3 


fid^vi = 


V- Vi 


d fi = -j—prR + RsM, 


i?2(v). 


where Ri and R 2 are functions that vanish faster than l/|v| for large |v|, while R 3 is 
a function that vanishes faster than l/|vp for large values of |v|. 


The proof of the theorem depends on the 

Lemma 3.2. Let /(v) be a function in §(K^), v G and h{v) := |v|/(v). The 
Fourier transform of the Green operator applied to h, namely the Fourier transform 

of 

(3.2) 9{h) := [ 

Jr3 |v- vi| 

is a function in the L^ space. Proof of Lemma 3.2. Since / is a Schwartz function, 
both h and Ah are in T^(]R^). For the Riemann-Lebesgue lemma their Fourier trans¬ 
form T{h) and T(Ah) are uniformly bounded. By standard properties of the Fourier 
transform, we have that 


(3.3) nH)K < 

for some positive constant C. Let’s recall now that the Green operator in Fourier 
space is a multiplicative operator namely T(S(h))(k) = 3’'(/i)(k)|k|“^, hence 


m9{h)){k)\ < 


C 

|kP(l + |kP) 


which means that |T(S(h))(k)| is bounded by an element of L^(R^). So, we have the 
thesis. □ 

We can now proceed with the proof of Theorem 3.1. Proof of Theorem 3.1 We 
will first prove 1) computing the limit 


lim |v| 

I v| —)-C30 


Jr3 |v-V l| 


d^vi 


lim 

I v| —>-oo 



|v| - |vi| 

|v-Vi| 


fid^vi + 




where we used the linearity of the Green operator in order to divide the desired limit 
in two parts. Let us discuss them separately. As shown in Lemma 3.2, the second 
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integral in the previous expression is the inverse Fourier transform of an fuirction, 
hence, by Riemann-Lebesgue lemma, its limit for large |v| vanishes. Since the absolute 
value of the function (|v| — |vi|)/|v — Vi| is bounded by 1 uniformly in v, the limit 
of the first integral can be obtained applying the dominated convergence theorem, 
namely taking the liirrit before computing the integral. The final result is 

lim |v| [ -——— rd^vi = n 
|v|->-oo Jr3 |v — Vi| 


and this proves the first assertion. 

Let’s proceed with 2). We observe that 

1^1 I" (v- Vl)'(v- ^ f (|v| - |Vl|) (v- vQ^V- VlV 




- vJ3 


V - Vl 


|v-Vi| 

|vi| (v - Vi)'(v - VrV 

|v-vi| |v-vip 


V- Vl 


fid^vi- 


fid^vi. 


For large values of |v|, the limit of the first term of the last sum can be obtained 
applying the dominated convergence theorem. The second term tends to zero as |v| 
is very large. In fact, as fi is positive. 


0 < 


|vi| (v - Vl)'(v - VlV 


V- Vl 


|v-Vi| 


fid^vi 


< 


|vi|/l ,3 


/r3 |v-V i| 


Vl- 


We can thus apply an argument similar to the one employed in the proof of 1) to 
show that it vanishes. Hence we have that 


lim 

|v|^-o< 


(v - Vl)^(v - Vl)-^ 




fid^vi = 


v'v' 


-n. 


from which follows the thesis. 

In order to prove 3) we observe that 


/r3 |v - Vl I 


id'^vi 


(|v| - |vi|) 

|v-Vi| 


djid^vi + 



djid^vi. 


Since / G §(R^), also its partial derivatives are in S(R^), i.e. dif G §(IR.^) for each 
f = 1, 2, 3, and we can thus proceed as in in the proof of 1). 0 

Applying the result of the previous theorem to (3.1) and recalling the discussion 
about the decaying properties oi Sb given at the beginning of this section, we have 
the following 

Proposition 3.3. The linearized Landau integral kernel Ql = Lb + Sb, com¬ 
puted over a background distribution f G which is of Schwartz type in the velocity 
domain, is such that 


SttAc^ / 71 71 \ 

QB{5f) = LB(t,x,v) + 5B(t,x,v) = + R. 

where R is function which vanishes faster than l/|v| for large values of |v|. 

We specialize now the asymptotic expansion provided by the previous proposition 
for backgrounds described by Maxwell distribution 


(3.4) 


fit, 


. n{mj5)'i _r 

,v) = — 


\/87r2 
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where /3 = (ki,T)~^, ki, is the Boltzmann constant and T the temperature of back¬ 
ground particles. In particular, we have that for large (3, i.e. in the limit of vanishing 
temperature (cold target), the asymptotic expression found above becomes exact, and 
so the following proposition holds. 

Proposition 3.4. The linearized Landau kernel Ql = LB{t,:x.^v) -|-S'b)!, x, v) 
defined as in (2.5) for perturbations Sf over a Maxwell distribution f as in (3.4) at 
inverse temperature /3 is such that 

OiW)(t,x,v) = + ipyap/) W i" «■ 


Proof. In order to prove this statement, it is sufficient to notice that, in the limit 
of large /3, the Maxwell distribution is proportional to a regularization of the Dirac 
delta centered in zero, where the factor of proportionality is noting but the plasma 
density. Hence, being Sb local in /, it vanishes for |v| 0. At the same time, 

analyzing Lb in (3.1) we have the thesis. □ 

In the previous proposition we have seen that the leading contribution in the 
linearized Landau kernel Ql for large values of |v|, is expressed by the following 
operator 


LB{t,x,v) 


SttAc^ 



n vV-^ 



We shall now specialize the form oi Lb for 6 f which are isotropic in v namely which 
are functions of |v|. 


(3.5) 


SfM = 5(|v|). 


For this kind of functions, by direct computation, it holds that 


(3.6) 


LB(t,x, v) 


Stt Ae^ n , SirAe^ n 

2 F] V2^ 2 F] 

2|vp 2 |vF 


3.2. Asymptotic form of the energy loss term and its relation with the 
Landau integral kernel. In the previous section we have analyzed the asymptotic 
form of the loss term (2.6) of the Landau equation (2.5) in the limit of high values of 
electron velocities, which is equivalent to the limit of small temperatures in the case of 
thermal background. Here we compare the asymptotic form of the Landau collision 
kernel obtained so far with a stochastic loss term added to a continuity equation 
(1.3) for the linearized perturbations 5f. Let us recall the equation for the linearized 
perturbation 5f 


^ + V • V^5f + - • Vv^/ 
at m 


jy^E) 

m 


Vv(5/-kS'(t,x,v), 


where (VxA) is the energy loss due to collisions. The loss term due to collisions 
in the previous equation is thus given by 

(3.7) Lc(^/)(Lx,v) := . v,5/. 

m 


Let us now briefly discuss the form of the stochastic force (Vx^l) due to Coulomb 
collision of the perturbed particles distributed as Sf with background particles in the 










ensemble /. We start recalling that the Coulomb collisional energy loss rate obtained 
in [36, 11, 28, 15] is 




vi\,a)ATfi(f<7d^vi 


where B is defined in terms of the Rutherford cross section as in (2.2) and AT is 
the variation of the kinetic energy of an accelerated electron due to a single collision, 
namely 


AT = - 


m|v — vi I 


-W • V — 


m|v — vi| 


-{cos9 — 1), 


where uj = a — (0,0,1) = (sin0cossin0sin</),cos0 — 1) is the deviation angle [11]. 
Performing the integral over all possible scattering angles and recalling that (^) = 
V • (VxT), the stochastic force results as being 


(VxT) = 


r f Vvi/i ,3 _ V 

m Ur3|v-vi| 2|v|2 


/i 

|v- Vi| 



From this expression it is easy to write the loss term Lq and, furthermore, its asymp¬ 
totic form can be obtained using Theorem 3.1. In particular the following proposition 
holds 

Proposition 3.5. Let f be any smooth function whieh is of Schwartz type in the 
velocity domain. The asymptotic expansion of Lc(t.,x,-v) defined in (3.7) for large 
|v| is such that 

Lc{t,yi,v) =Tc(i,x,v) +R = Vvd/+ T 

where, for large values o/ |v|, R is a function that vanishes faster than l/|v|^. 
When the background is thermal, namely when / is a Maxwell distribution as in 
(3.4), proceeding as in the proof of proposition 3.4 we have that the asymptotic form 
obtained in the previous proposition implies that the limit of vanishing background 
temperature is such that 

lim Tc = Tc = • Vv(5/, v ^ 0. 

/3->oo 2|v|'^ 

Finally, by direct inspection, we notice that whenever Sf is an isotropic perturba¬ 
tion as in (3.5), the asymptotic form of the stochastic loss term of the continuity 
equation Tc(t,x, v) matches the asymptotic form of the Boltzmann collision kernel 
Lsit, X, v) given in (3.6). This observation suggests that the correct form of the mod¬ 
ified continuity equation is (1.3) and not (1.4) as sometimes assumed in the literature. 
In particular, the stochastic force due to Coulomb collisions needs to be taken into 
account adding to the continuity equation a contribution of the form 

(3.8) Lc{Sf)(t, X, v) = ■ y^6f. 

m 

4. Resolution of the continuity equation with hyperbolic methods. The 

aim of this section is to write the continuity equation and its solution in a form suit¬ 
able to be compared with data obtained by satellite imaging-spectrometer observing 
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solar flares. At this point it is important to recall that these instruments measure 
photon fluxes resolved in space, time and energy However, since photons are produced 
by accelerated electrons due to Bremsstrahlung, it is in principle possible to directly 
obtain information about electron ensemble from the measured photon fluxes. This 
has been accomplished by Piana et al. in [32] where they employed refined regular¬ 
ization techniques to invert the Bremsstrahlung cross section. The outcomes of the 
analysis presented in [32] are maps of electron fluxes F(t^ x, A), which are noting but 
images resolved in space and (kinetic) energy. 

It is at this point interesting to compare theoretical predictions for electron densi¬ 
ties /, following from certain kinetic equations like (1.3), with observations of electron 
fluxes F obtained from satellite measurements. 

However, since electron fluxes are integrated quantities over all velocity directions, 
this comparison is not completely straightforward. 

For this reason, starting from the precise relation among F and /, we transform 
(1.3) into an equation valid for the electron flux F{t, x, E). Afterwards we shall solve 
this equation using some hyperbolic techniques to obtain a solution comparable with 
measurements. 

Let us recall that, at fixed time, F describes the flux density over the domain 

X 1R+ parametrized by (x, A), while / describes the particle density over R® 
parametrized by (x, v). Since E is the kinetic energy of a particle of velocity v, 
the definition of F in terms of / is 

(4.1) F(t, X, if) ^ / /(t, X, v)dS2 , E = 

m Js,2 2 

where the integration is over all possible velocity directions and hence is the 
standard measure of the two-dimensional sphere Similarly, having in mind (1.3), 
the same averaging needs to be applied on the source term. Then, we can define the 
injection flux S related to the source S to be 

(4.2) S(t,x,F;) := / S'(t,x,v)d§2. 

m i§2 

Although the information about the velocity directions is lost in F, under some par¬ 
ticular assumptions we are going to discuss, we might be able to extract an equation 
for electron fluxes F from (1.3). 

First of all, in the limit of small background temperature, the leading contribution 
in the electron distribution /(t,x, v) is given by accelerated electrons only. For this 
reason, the thermal background might be discarded and thus / obeys the continuity 
equation equipped with the stochastic force (1.3), 

df 1 

(4.3) ^+vVx/ =- (VxL;) • Vv/ + S{t, X, v) 

ot m 

where 5'(f,x, v) is some external source and (V^E) is the energy loss rate. 

Magnetic fields are particularly strong during solar flares. We thus assume that 
electrons move freely along the magnetic field lines discarding the rotational motion 
in perpendicular planes. At the same time we expect that the electron distribution / 
does not vary much passing to nearby magnetic field lines. Under this assumption it is 
convenient to parametrize the space with an appropriate coordinate system (s,a;,y), 
and hence s is the arc-length of a magnetic field line 7 while djdx and d/dy are 


10 


orthogonal to 7 . Thanks to this choice f{t, x, v) = f{t, s, v), namely the distribution / 
depends on space only through s while it is constant under changes of x, y coordinates. 
In order to obtain an equation for the flux density F, we assume that both / and 
the source S in (4.1) and in (4.2) depend on the velocities only through |v|, hence 
/ = |v|). 

We proceed writing the accelerated electron ensemble as the sum of two parts 
/ = /+ + /_, where particles in /+ move towards positive coordinates s while those 
in /_ towards negative ones. In other words, /+ (resp. /_) vanishes when Vg, the 
component of the velocity v along 7 , is negative (resp. positive). Furthermore, (4.3) 
holds separately for /+ and /_, hence 

(4.4) =_l(V,F;).Vv/± + ^(t,x,|v|). 

ot os m 


To obtain an equation for the flux densities F± corresponding to /±, we proceed by 
computing the integral of the previous two equations on the velocity directions as in 
(4.1) and in (4.2). Thanks to the isotropic hypothesis for / and S, the right hand side 
of (4.4) does not depend on the velocity directions, while the left hand side depends 
on the directions only through Vg. Since the average value of Vg over the velocity 
directions, is |v|/2 (resp. —|v|/2) for particles in /+ (reps. /_), the integrals gives 
(4.5) 


d F±{t,s,E) 
dt E 


± 


d F±{t,s,E) , 

/dE\ 

^ 1 

^F±(t,s,F)\ 

V 2 m ds E ' 

\ ds / 

^ V m dE ' 

k E J 


F,±{t,s,E) 

E 


where T,j^{t,s,E) (resp. S-) are obtained by an integration of S{t,s,v) as in (4.1) 
over velocity directions corresponding to positive (resp. negative) Vg. Notice that, 
since S is assumed to be isotropic in v, we have that 11+ = I]_. 

We shall now look for static solutions of the previous equation. Assuming = 
E_ = E(s, E) to be constant in time, (4.5) reads 


^ d F±is,E) ^/dE\ d E±{s,E) /^ E(s,A) 

ds E \ds / dE E \ E E 


We shall now simplify the problem by means of the following substitutions^: 


(4.7) A{s,E) = 


E E 


x{E) X + f ^ / dE\ ’ 
JE 2{-^) 


F± 


Then (4.6) take the form 

(4.8) ± —4>±-—$±=A. 

Since a constant solution is not physically admissible neither for F^ nor for FL, we 
can apply on both sides of (4.8) the operators 


ds dx 

without loosing any information. The mixed terms disappear and thus both equations 
for <i>+ and <i>_ can be combined to give 


(4.9) 


□ $ = 2—A 
dx 


^The definition of x is obtained imposing ^ = 2 



and supposing x{E) = x 
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where $ = $+ + $_ = F/E is the sum of forward and backward moving electrons 
and □ is the D’Alembert operator in two dimensions. Equation (4.9) is nothing but 
the wave equation in two dimensions with a source term, where x plays the role of 
time and s of space. 

To find the solution of this equation we have to provide physically motivated 
initial conditions. Since electrons cannot gain energy during their motion (injected 
electrons can only loose energy during their path), we shall look for solutions which 
satisfy the following initial conditions 

(4.10) lim $(a;, s ± x) = 0 , lim dx^{x, s ± a;) = 0 , Vs . 

x—>-oo x—^oo 

The solution can now be obtained by means of the convolution (in the x, s space) of the 
advanced (with respect to x) fundamental solution (Green function) 6- = 5{x |s|) 

of the D’Alembert equation with the source term 2A, namely 

^ = 5-*2A. 


The absence of the derivative dx in front of A in the previous expression, contrary 
to the case of (4.9), is due to the fact that □(!_ = dx5. Actually, in two dimensions 
such a derivative is necessary in order to avoid infrared divergences. Furthermore, 
F obtained in that way is the unique solution obeying the desired initial conditions 

(4.10) . All other solutions of the wave equation with the chosen source do not vanish 
for large energies. 

We can now perform the convolution using (4.7), namely 

(4.11) F(s, E) = 2E J (5 (a; - xo + |s - Sol) ^ ^ ^ ^ ’ 

where x{E) is defined as in (4.7) and E{xo) is the inverse function of x{E) evaluated 
in xq. So, (4.11) gives 


(4.12) F{s,E)=2E 


2m 


S(so, E{x + |s - Sol)) 


i?(x + |s-so|) E(x+|s-so|) 


dsn 


Notice that the solution written in the form (4.12), is in accordance with the one 
presented in [20] under the assumption of uniform injection in a limited region of the 
flare. Furthermore, this expression can be directly compared with results provided by 
experiments, as in [39, 12]. 


More precisely, as in the case of cold target [15], the energy loss rate is obtained 
from 



—27re^An 



letting AT = 27r A, we have 

^2 

x(E) = - and hence E(x) = 2\/Knx, 

AKn 

and (4.12) takes then the form 


(4.13) F{s,E)=2E 


2m 


(F2 +4nA:|s 


_ S(so,(F^+4nA-|s-so|)^^^) 

sol)^/" {E^+AnK\s-so\f/^ 
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5. Comparison with experiments. In this section we validate the solution of 
the continuity equation obtained so far for accelerated electron fluxes with the results 
of a real measurement occurred during a solar flare. First of all, we compare the 
outcome of our analysis with the continuity equation for the electron flux described 
by Emslie, Barrett and Brown [16]. Following the notation of the previous sections 
this equation reads 


(5.1) 


dF{s,E) 

ds 


A 



F{s,E) 


^{s,E). 


Applying the same procedure used to solve (4.5), we obtain a solution for (5.1), in 
the form of 

(5.2, f (., E) = 2 (f)" / {f(Ei.'))) 

where x' = x \s — sqI and, in contrast to (4.7), x is defined as 

(5.3) x{E)=x+ f j^dE'. 

•>E (- 37 ) 

From a theoretical point of view, the expressions of E{s,E) given by (4.12) and (5.2) 
differ by two aspects. Firstly, the energy loss rate enters (4.12) and (5.2) in two dif¬ 
ferent ways. The resulting electron fluxes are compared in Figure 1 where we notice a 
significant difference at lower energies, e.g. 1-2 keV. Secondly, (4.5) is deduced from 

(4.3) performing the integral over all particle directions. In this way, collisions along 
every direction are taken into account. On the contrary, in the derivation of (5.1) only 
interactions along the magnetic field lines are considered. These different assumptions 
have an explicit effect on the expression of x{E), used in (4.7) and in (5.3), making 
them differ by a factor two. This modification is evident also at higher energies, in 
particular at 10 — 100 keV, the energy range in which measurements can be detected 
by imaging-spectrometers. 


We check now that the electron flux distributions described in (4.12) and in (5.2) are 
both in accordance with the available experimental data. To accomplish this task 
we recall that images of electron fluxes deduced by the NASA Reuven Ramaty High 
Energy Solar Spectroscopic Imager (RHESSl) spectrometer [26] can be analyzed. We 
shall thus validate both solutions (4.12) and (5.2) with the observations of a single 
event. Furthermore, regarding the form of the source considered in the theoretical 
solutions, we assume that the injected electrons are constant along a certain region 
of the flare with length Lq be. 

(5.4) s(»,B) = k,^(|-) 

where 0 is the Heaviside function, hs is a constant which represents the intensity of 
the injection, Eq is the minimum injection energy and Sm represents the origin of the 
arc-length system of coordinates. Moreover, the hot target model described in [36] is 
used as energy loss rate, i.e.. 



(5.5) 
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Comporison between continuity solutions 
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Fig. 1. In red and black the profiles of g{E) obtained respectively with F{s,E) measured as in 
(4.12) and in (5.2) in the electron energy range between 1,75 keV and 2 keV. 

where n represents the background electron density, kT is the plasma temperature 
measured in keV and erf is the ordinary error function. As experimental data, we 
consider the RHESSI acquisitions of a flare, which occurred on the 15th of April 2002 
in the time range between 00:00 and 00:05. 

More in details, we have selected an electron energy range between 18 keV and 
34 keV. We have separated this interval into eight sub-intervals of amplitude of 2 
keV and we have recovered the correspondent electron flux. In particular, using the 
method described in [32], we have obtained eight electron flux maps similar to the 
one shown in the first column of Figure 2. 

As discussed in the theoretical analysis presented in §4, we have to restrict the 
attention to the motion of accelerated electrons along a single fixed magnetic field 
line. The magnetic field line, or better the electron path, we have considered has 
thus been selected as the locus of the brightest points of each image row, following 
the procedure described in [39]. The value of each pixel along this path represents 
the measured electron flux F{s, E). The chosen pixels and the relative electron flux 
values are shown in the second and the third columns of Figure 2 respectively. 

Then, for each energy channel, we have computed the mean value of the electron 
flux along the path, obtaining an averaged electron flux g{E) of the form 



(5.6) 


The measured values of g{E) are shown in the first row of Figure 3. 

We have thus deduced this same quantity using the theoretical models described 
in (4.12) and (5.2) for the electron flux, provided a source term and an energy loss 
rate as described in (5.4) and (5.5). Hence, the measured values of g{E) have been 
fitted against theoretical ones using a simulated annealing algorithm [22] with the 
mean squared error percentage as an energy function. The best values for hg and 
S have been obtained, while the target electrons density, the plasma temperature. 
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the minimum injection energy and the injection amplitude have been fixed as n = 
10^^ cm“^, fcT = 1, 75keV, Eq = 18keV and Lq = SOarcsec respectively. 



Fig. 2. In the first column a reconstruction of an electron flux map obtained by the algorithm 
described in [32] starting from RHESSI data is provided. The accelerated electron path selected 
following the method described in [39] is highlighted in white in the second column. The associated 
pixel values (i.e. the electron flux amounts) are shown in the third column. 

The values obtained are shown in table 1 while in Figure 3 the experimental pro- 
file and the fitted solutions are compared. It’s interesting to notice that both (5.2) 
and (4.12) fit very well the measured data, as both models give a really low value 
for the mean squared error percentage. Furthermore, the parameters obtained are 
physically consistent and, actually, they are all similar to those found in [20, 12, 39]. 
On the other hand, although the fitted curves are almost indistinguishable, we stress 
that parameters obtained using the two models are slightly different. In principle it is 
thus possible to distinguish between (4.12) and (5.2). These difference should become 
more evident when more refined measurements will be available. 

Summarizing, we have compared two different models for the theoretical descrip¬ 
tion of the electron flux, described by (5.2) and (4.12) respectively. It is important 
to point out that while the former model has been deduced under an hypothesis of 
flux conservation, the latter has been developed starting from basic principles and it 
is still in accordance with experimental data. 


hs [cm ^ keV 

<5 

MSEP 

1.39 X 10’' 

6.81 X 10® 

8.12 

7.88 

0.0003 

0.0006 


Table 1 


Plasma parameters obtained with (4.12) and (5.2) respectively in the first and in the second 
row: in the first column the obtained injection intensity, in the second the amplitude of the injection 
region, in the third the injection spectral index and in the fourth the resulting mean sguared error 
percentage. 
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15th April 2002 00:00-00:05 



15th April 2002 00:00-00:05 



Fig. 3. 2002 April 15 event: g{E) as a function of electron energy. First row: Empirical values 
of the averaged electron flux. Second row: Empirical values fitted against theoretical models, in red 
the profile obtained with (4.12) and in black the profile obtained with (5.2). 


6. Comments and open issues. In this paper we showed that in some asymp¬ 
totic limits the linearized kinetic Boltzmann equation for Coulomb collisions reduces 
to a continuity equation equipped with an external stochastic force. Specihcally, un¬ 
der the assumption of grazing collisions (namely when the Coulomb logarithm is very 
large), we showed that, in the limit of high values for the electron velocities, the Lan¬ 
dau equation reduces to the continuity equation, provided a stochastic loss term of 
the form 


(V^x) 

m 


•V/, 


is taken into account. The resulting equation has the form 




Vv/ + 


{V^E) 


^.f = s. 
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We point out that the stochastic contribution written above, is slightly different than 
the usual loss term considered in some literature about electron motion in solar flares, 
especially at lower values for the electron energy E. Since measurements in these 
domain start becoming very precise this difference should be appreciated in the close 
future. For this reason, in the second part of the paper, we have rewritten the con¬ 
tinuity equation in an equivalent form ready to be compared with data collected 
by spectrometer measuring X-ray photons emitted during solar flares. Indeed, the 
continuity equation has been reduced to an hyperbolic equation allowing a simple 
closed form solution. Finally, the obtained solution is written in a form which can 
be compared with experimental data collected by satellite hard X-ray spectrometer 
analyzing the sun, see e.g. (4.12). In the last section we have shown that the changes 
suggested for the description of the motion of accelerated electrons are still in accor¬ 
dance with data recorded by the RHESSI spectrometer. More generally, it has been 
shown in [16, 39, 12, 20] that the continuity equation can be used as a tool for the 
validation of theoretical models for the injection mechanisms and for the accelerated 
electrons energy loss using data provided by different imaging instruments. For this 
reason, having a more manageable equations is important, as it allows a more accurate 
analysis to choose among different acceleration models. 
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